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Abstract 

A unified view is given to recent developments about a systematic 
method of constructing rational mappings as ergodic transformations 
; with non-uniform invariant measures on the unit interval / = [0,1]. 

| All of the rational ergodic mappings of / with explicit non-uniform 

■^J- ■ invariant densities can be obtained by addition theorems of elliptic 

functions. It is shown here that the class of the rational ergodic map- 
\ pings I — > I are essentially same as the permutable rational functions 

obtained by J. F. Ritt. 
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1 Introduction 



Qualitative characterizing discrete-time dynamical systems have been recently 
gaining an attention. In continuous-time dynamical systems given by ordinary 
differential equations (ODE), the notion of integrability can serve as the sharp 
characteristic for general ODEs such that integrable ODEs are solvable in the 
sense that exact solutions are analytically obtained and non-integrable ODEs 
such as the three-body problems are not solvable and they show, in general, 
chaotic behavior. In contrast with ODE cases, the notion of integrability itself 
is still vague in discrete-time systems, that is, it is not so sharp as the case of 
ODEs. Quite recently, Hietarinta and Viallet show examples of discrete-time 
dynamical systems which pass the singularity confinement test (discrete-time 
version of Painleve test), but which nevertheless show chaotic behavior [|IJ. 
Thus, an issue is now focused on searching a characteristic distinguishing 
chaotic discrete-time systems from regular discrete-systems (which have been 
thought as integrable systems. The purpose of the present paper is to present 
a class of discrete-time dynamical systems which have exact solutions, but 
which nevertheless are shown to be ergodic (thus, chaotic). Such dynamical 
systems (we call exactly solvable chaos) can be constructed in a systematic 
way, by utilizing addition theorems of elliptic functions. The main claim here 
is that solvable discrete-time dynamical systems with chaotic properties do 
not belong to exceptional classes but they are as ubiquitous as multiplication 
formulas of elliptic functions, which have interesting implications to the issue, 
on the issue. 



2 Construction 

Let us consider an elliptic function s(x) whose inverse function is defined by 
the following elliptic integral: 

i. , r x du r x2 dv 

Jo Vao + a2U 2 + a 4 u* + a 6 « 6 Jo 2^v(a + a 2 v + a 4 v 2 + a 6 v 3 ) ' 

(1) 

where < x < 1 and a , ■ ■ - a^ are real-valued constants, which satisfy the 
following condition: 

a + a 2 u 2 + azfU 4 + a 6 w 6 > for < u < 1 , , 

do + «2 + «4 + 0-6 = 0, a>o > 0. 

This elliptic function can be also considered as the solution q(t) of a Hamil- 
tonian system the motion with a Hamiltonian 

H = \? + V(q) = 0, (3) 
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where the potential function V(q) is an even polynomial function of the form 
V(q) = —(cio + a 2<l 2 + tt^q 4 + a 6 g 6 ) satisfying the conditions 

V(0) = 0, V(l) = 0, V(q) < for < q < 1. (4) 

Thus on this condition about a , • • • , in Eq. (0), this elliptic function s(x) 
has a real period 4K, where K is defined by the elliptic integral 

K= f 1 du (5) 

Jo \/ r ao~+ CL2U 2 + a^u 4 + a 6 w 6 

We have the relations 

< s{x) < 1 for < x < AT, (6) 

and 

s(0) = 0, s{K) = 1, s(2if) = 0, s(3K) = -1, s(4/f) = 0. (7) 
Since the equality 

S_1 (-^) = / a. 2^ 4 4 x 6 6 = -^^^ S 2/ (8) 

holds, then we have 

s(y) = -x, s(-y) = x (9) 

which means that s 2 (y) is an even function as s 2 (y) = s 2 (—y) and has a real 
period 2K. Furthermore, since any elliptic function has an algebraic addition 
theorem ||, there exists a polynomial function F in three variables such that 
the relation 

F(s 2 (x! + x 2 ), s 2 (xx), s 2 (x 2 )) = (10) 

holds. As a special case of the addition theorems, there exists a polynomial 
function G in two variables such that 

G(s 2 (px),s 2 (x)) = 0, (11) 

where p is a positive integer greater than the unity. This means that we have 
an algebraic mapping as 

s 2 (px) = f(s 2 (x)), (12) 
where / is an algebraic function. If we set 

X n+1 = s 2 (px), X n = s 2 (x), (13) 

we have a discrete-time dynamical system 

X n+1 = f(X n ). (14) 
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on the unit interval I = [0, 1]. It is easy to check that there exist p pre-images 

0<x m = s 2 [ < 1, 0<m<p-l 15 

V 

satisfying y = f(x m ) for < y < 1. This dynamical system X n+1 = f(X n ) 
has the following remarkable property: 

Theorem 1 A dynamical system X n+ \ = f(X n ) defined as the addition the- 
orem ftT^) of an elliptic function s 2 (x) is ergodic with respect to an invariant 
measure n(dx) which is absolutely continuous with respect to the Lebesgue 
measure and their density function p(x) is given by the formula 

p( x ) = r=, 1 ( 16 ) 

2KJx(ao + a2X + a^x 1 + a§x 6 ) 

(Proof of Theorem [I]) 

Let us consider the diffeomorphisms of / = [0, 1] into itself given by 

< (j){x) = — s~ x {^x) < 1 for < x < 1. (17) 
K 

Using the relations 

s\K-p9) = f[s\K9)} for 0e[O,i], 
s 2 [K{2-p6)]=f[ S \K6)] for 9 G [i,2], 

s^[K(-(-l)H+ 1 -^y + (-iyp9)] = f[s\Ke)] for £P € [j, 

s "[K((-inp-l) + 1 -±^ + (-iy- 1 p9)] = f[s\Ke)} for Jg[^,1] 

(18) 

by the fact that s 2 (x) is an even function and it has a real period 2K, we can 
derive the piecewise-linear map f(x) = <p ° / ° <p~ 1 {x) on I = [0, 1] as 

f(x) = px for x G [0, -] 

f( x ) = -(-i)H+ 1 -^ + (-iypx xe[y-f\ 
fix) = -i-iy-Kp - 1) + + (-i) p - V x g [2=1, i] 



(19) 



It is noted here that when p = 2, /(x) is the tent map, given by 

f(x) =2x x G [0, §] 
/(x)=2-2x xG[|,l]. 



(20) 



Clearly, the map f(x) in Eq. fllTf ) is ergodic with respect to the Lebesgue 
measure on I = [0, 1]. Thus, / is also ergodic with respect to the measure 

Ai(dx) = ^dx = ±- dX (21) 



The measure fi(dx) is absolutely continuous with respect to the Lebesgue mea- 
sure, by which we can define the algebraic density function p(x) in Eq.flTEl). 
(End of proof) 

A simple corollary of the theorem is the following: According to the Pesin 
identity which can be applied to the dynamical systems with absolutely con- 
tinuous measure ||, the Lyapunov characteristic exponent A(/x) for the mea- 
sure of the map A n+1 = f(X n ) is equivalent to the Kolmogorov-Sinai entropy 
h(p) = logp, namely: 

rl df 

AO) = J o log |^| • p{x)dx = h(n) = log p. (22) 

In this sense, discrete-time dynamical systems X n+ \ = f(X n ) constitute 
a typical class of chaotic dynamical systems with a special property that 
their invariant density functions are explicitly given. Hence, we call this class 
of chaotic dynamical systems exactly solvable chaos [|J. We note here that 
dynamical systems of exactly solvable chaos can be good pseudo random- 
number generators |J, especially for Monte Carlo simulations from the exact 
ergodicity 

N-l j ! 

i im TfQM = / Q{x)p{x)dx. (23) 

We remark here that we can generate infinitely concrete chaotic dynami- 
cal systems with the unique invariant measure (|T6| ) from addition theorems 
s 2 (px) = f[s 2 (x)} of elliptic functions s(x) for p = 2, 3, • • •. 



3 Examples of exactly solvable chaos 

Historically, the most simple example of exactly solvable chaos is the logistic 
map Y = 4X(1 — X) = fo(X) on I = [0, 1] given by Ulam and von Neumann 
in the late 1940's§. 

Ulam and von Neumann show that A n+1 = 4X n (l — X n ) is ergodic with 
respect to an invariant probability measure a(dx) = — , dx The Ulam- von 

-KtJ x{\—x) 

Neumann map can be seen as a special case of exactly solvable chaos of the 
above theorem 1 when we consider s(x) defined by 

s^ix) = sin" 1 ^) = [ J U 2 . (24) 



o VT 



u 
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This corresponds to the case that 

do = 1, 0,2 = — 1, CI4 = Gt6 = (25) 

in Eq. (0). Thus, we have the relation s(x) = sin (a;) which admits the 
duplication formula sin 2 (2x) = 4sin 2 (a;)(l — sin 2 (a;)) = /(sin 2 (x)) (p = 2 in 
Eq.(^)) gives the logistic map. Clearly, in this case, K — |-. Thus, using 
Theorem [I], we have calculate the density function p{x) as 

= / * ■ (26) 
nyx{l — x) 

When p = 3, we have i/ie cu&ic map as X n+ i = X(3 — AX) 2 from the tripli- 
cation formula 

sin 2 (3x) = /[sin 2 (x)] = sin 2 (a;) [3-4 sin 2 (a;)] 2 (27) 

The cubic map can be regarded as a special class of Chebyshev maps obtained 
by using more general addition formulas as sin 2 (px) = /(sin 2 (px))|^]. 

The alternative attempts to generalize the Ulam-von Neumann map within 
a set of rational functions was made by Katsura and Fukuda in 1985 M. The 
Katsura-Fukuda map is given by 



4X(1-X)(1-IX) 
(1 - /X 2 



Y = ^LZ^Zl^I s fl{x) {2t 



for < I < 1. Clearly, the Ulam-von Neumann map can be regarded as a 
special case of Katsura-Fukuda maps. In this case, the corresponding elliptic 
function s(x) is Jacobi sn function whose inverse function is defined by 



du 

-M 2 )(l -lu 2 ) 



s - l (x) = sn- l (x;Vl) = / . , (29) 



where Vl corresponds to the modulus of Jacobi elliptic functions. This cor- 
responds to the case that 

ao = 1, a,2 = —I + 1, , Q>4, = I, «6 — (30) 

in Eq. ([!]). Katsura and Fukuda show M that the map ( PB| ) has exact 
solutions X n = sn 2 (2 n 8; \fl). Using the idea of theorem 1, the author shows 
H that the Katsura-Fukuda maps (p8|) are also ergodic with respect to an 
invariant measures which can be written explicitly as 

dx 

fi(dx) = p(x)dx = , , (31) 

2K(l)Jx{l-x)(l-lx) 



du 



where K(l) is the elliptic integral of the first kind given by K(l) = — , . 

yf (1 — U 2 )(l — lu 2 ) 

This can be easily checked using the duplication formula || of the Jacobi sn 
elliptic function 



sn(2 ■ Vl) = 2SU{U] ~ SU2{U] ^ m " lSn2{U] (32) 

Recently, the more general classes of exactly solvable chaos are derived 
from elliptic functions s(x) whose inverse functions are defined by 

s-\x) = [ X dU (33) 

Jo ^(1 -u 2 )(l -lu 2 ){l -mu 2 ) 

where the parameters I and m are arbitrary real numbers satisfying the con- 
dition — oo < m < I < 1 Hj. This corresponds to the case that 

a = 1, a 2 = —(I + m + 1), a 4 = lm + I + m, a 6 = — Zm (34) 

in Eq. ([TJ). The associated mapping generated by the duplication theorem 
s 2 (2x) = fi 2 m(s 2 (x)) is given by the following rational transformations 

(2) 4I(l-I)(l-II)(l-mI) 
hmW- 1 + AX 2 + BX 3 + CX 4 Gi ' {6b) 

where A = —2(l + m+lm), B = 8lm,C = I 2 + m 2 -2lm-2l 2 m-2lm 2 + l 2 m 2 , 
and X € Clearly, this mapping is a generalized version of Ulam-von 

Neumann map and Katsuda-Fukuda map. This two-parameter family of the 
dynamical systems x n+ i = /j^v^n) (|35Q are also ergodic with respect to an 
invariant measures given by 

n(dx) = p(x)dx = , (36) 

2K(l, m)\Jx(l — x)(l — lx)(l — mx) 

where K is given by the integrals 

K(l, m) = C = (37) 

J o ^(1 -u 2 ){l - lu 2 ){l - mu 2 ) 

We can check this fact by directly computing the duplication formula of s(x) 
whose inverse function is defined in Eq. fl3"3"|). Let us represent s(x) in terms 
of the Weierstrass elliptic functions 

The Weierstrass elliptic function p(u) of u G C is defined by 

1 ' 1 1 

P{U) = + 5 2M-2ku 2 y ~ (2juj 1 + 2kuj 2 ) 2h (38) 
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where the symbol J2' means that the summation is made over all combinations 
of integers j and k, except for the combination j = k = 0, and 2uj\ and 2u 2 
are periods of the function The Weierstrass elliptic function p(u) 

satisfies the differential equation 

(^) 2 = V(x)-^p(x)-^ (39) 

with the invariants 

(40) 

Let ei, e2 and e3 be the roots of 4z 3 — — (73 = 0; that is, 

ei + e 2 + e 3 = 0, eie 2 + e 2 e 3 + e 3 ei = eie 2 e 3 = ^. (41) 

The discriminant A of the function p(u) is defined by A = g\ — 21g\. If 
A > 0, all roots ei,e 2 and e 3 of the equation 4z 3 — g 2 z — g 3 = are reaZ. 
Thus, we can assume that t\ > e 2 > e 3 . In this case, to>i and o; 2 are given by 
the formula 

r°° dz r e z dz 

u>i — -. , lu 2 = i / ■> ■ (42) 

J ei ^/Az i - g 2 z - g 3 J-00 ^g 3 + g 2 z - Az 6 

Using a rational transformation of a variable as v = ^+21-^^ — + 1, we 

can rewrite s" 1 (x) as 

2l + 2m-3lm-l . m) , 

i 4(l-l+l 2 -m+m 2 -Im) 4:(2-l-m)(2l-m-l)(2m-l-l) llr , n 

where # 2 = — ^ ^,#3 = — — ^ — -. We note here 

that 4y 3 — g 2 y — g 3 can be factored as 

o i( , 2-l-m.. 21 — m — l w 2m — Z — 1 
4y 3 - p 2 y - (? 3 = 4(y ){y ){y ). (44) 

Thus, we set e x = > e 2 = > e 3 = Thus, s(x) can be 

written explicitly in terms of the Weierstrass elliptic function as 

(l-Z)(l-m) sn 2 (Vl -mx; 
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s 2 (x) = 1 



pfa -X)- 2'+2m-3to-l 1 _ m + msn 2( v ^T^ x; ' 

145) 



S 



The function s 2 (x) also has the same periods 2u>i and 2lv 2 of the Weierstrass 
elliptic function p(x) and the real period 2u% is given by the formula 

2K( 

2K(l, m) = 2cux = . Kl - m \ (46) 
Vl — m 

Using the addition formula of Jacobi sn function, we finally obtain the explicit 
duplication formula of s(x) as 

As 2 (x)(l-s 2 (x))(l-ls 2 (x))(l-ms 2 (x)) 
{ ' l + As A (x) + Bs e (x) + Cs 8 (x) ' 1 ' 

where A = —2(l + m + lm), B = 8lm and C = l 2 + m 2 — 2lm — 2l 2 m — 21m 2 + 
l 2 m 2 . Thus, we obtain the two-parameter family of generalized Ulam-von 
Neumann maps X n+ i = //^ (X n )(|35|). 

(3) 

In the same way, we can construct generalized cubic maps from the trip- 
lication formula s 2 (3x) = fi 3 m(s 2 (x)) as 

X(-3 + 4X+_E<=i AiX^ 2 

^9 

where Ai, ■ ■ • , A 4 and B 2 , ■ ■ ■ ,Bg are given by 



A 1 = 4(1 + m),A 2 = -6(1 + m + lm),A 3 = 12lm, 
A 4 = l 2 + m 2 - 2lm - 2l 2 m - 2lm 2 + l 2 m 2 , 
B 2 = -12(1 + m + lm), 

B 3 = 8(1 + m + I 2 + m 2 + l 2 m + lm 2 + 15/m), 

B 4 = 6(5Z 2 + 5m 2 - 26/m - 2Ql 2 m - 26/m 2 + 5Z 2 m 2 ), 

B 5 = 24(-2/ 2 - 2m 2 - 2Z 3 - 2m 3 + Aim + 7l 2 m + 7lm 2 ) 

+2A(Al 3 m + 4/m 3 + 7Z 2 m 2 - 2/ 3 m 2 - 2/ 2 m 3 ), 

B 6 = A(Al 2 + Am 2 + Al A + 4m 4 + 17Z 3 + 17m 3 - Sim) 

+A(-lll 2 m - 17/m 2 - 17l 3 m - 17/m 3 - 8/ 4 m - 8/m 4 ) 

+4(4/ 2 m 4 + 4/ 4 m 2 - 17/ 3 m 2 - 17/ 2 m 3 + 17/ 3 m 3 - 54/ 2 m 2 ), 

B 7 = 2A(—l 3 —m 3 — I 4 — m A + l 2 m + lm 2 — l 3 m — lm 3 ) 

+2A(l A m + lm A + Al 2 m 2 + Al 3 m 2 + Al 2 m 3 ) 

+2A(l A m 2 + l 2 m A - l 3 m 3 - l A m 3 - l 3 m A ), 

B 8 = 3(3/ 4 + 3m 4 + 4/ 3 m + 4/m 3 + 4/ 4 m + 4/m 4 - 14/ 2 m 2 ) 

+3(-4Z 3 m 2 - 4/ 2 m 3 - 4/ 3 m 3 - 14/ 4 m 2 - 14/ 2 m 4 + 4/ 4 m 3 + 4/ 3 m 4 + 3/ 4 m 4 ), 

B 9 = 8(-l A m - lm A + l 3 m 2 + l 2 m 3 + l A m 2 + l 2 m A - 2l 3 m 3 + l A m 3 + l 3 m A - I 4 

The generalized cubic map fl^ has the same invariant measure ( pq) as 
the one of the generalized Ulam-von Neumann map. Furthermore, since 

s 2 (pmx) = flZ\s 2 iP2x)) = /<£> o fl%\s 2 (x)) 

= /&V(Pi*)) = fS ° fi%\s 2 (x)) = flZ P2 \s 2 (x)), 
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we have the following commutative relations 

JLm l A J - Jim ° Jim l A J - JLm ° Jim l A J> 



(49) 



where p% and p 2 are positive integers and f[)l(X) = X. Thus, from the func- 
tional relations (^), all generalized models of exactly solvable chaos given by 
Y = ftm(X) f° r & e can be constructed by generalized models of exactly 
solvable chaos given by Y = f^{X) with p being primes. For examples, we 
can compute Y = fi^l(X) as follows: 

fi%( X ) = fi% ° flfn( X ) = fll ° fl 2 m ( X ) 
fl!l( X ) = fll ° f^l ° fl 2 l( X ) 

/£(*) = /$ ° fi%m 

tr~ Pk \ x ) = /&> ° fti o • • • ffc\x). (so) 

All of the above examples of exactly solvable chaos are rational mappings on 
the unit interval / and they can be summarized in the following table. 



Explicit Invariant Densities of Rational Ergodic Mappings 



Rational Mappings 



Algebraic Densities 



Ulam-von Neumann map 
Chebyshev maps 
Katsura-Fukuda map 
Generalized Ulam-von Neumann maps 
Generalized Chebyshev maps 



p(x) 
p(x) 
p(x) 
p{x) 
p{x) 



w\J x(1—x) 
1 I 

1T-\J x(l — x) 



K(l)^x{l-x){l-lx) ° 
ftT(«,m) v /x(l-x-)(l-/a;)(l -rai 

1 



K 



(l,m)\/ x(l—x)(l—lx)(l—rnx) 



4 Ritt and Weierstrass's theorems 

The commutativity 

fp ° fq( X ) = fq° f P ( X ) = f P q( X ), ( 51 ) 

for rational ergodic mappings {f p } comes from the fact that their mappings 
are derived from the multiplication formulas of a single elliptic function s(x) 
as s(pu) = f p [s(u)]. Thus, they have exact solutions of the form X n = s(p n 9). 
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Thus, the commutativity is the key for one-dimensional ergodic dynamical 
systems to have exact solutions. 

With regard to the class of rational permutable functions in this way, Ritt 
showed in 1923||11|| that rational functions constructed from the multiplica- 
tion formulas of elliptic functions such as the above examples {f p } are all of 
rational commuting self-maps of the projected line except trivial cases. The 
trivial cases indicate the cases that rational permutable maps corresponds to 
the multiplication formulas of rational functions of e z , cos(z). Consequently, 
permutable rational maps which are shown to be ergodic in the previous sec- 
tions are essentially all of rational permutable ergodic mappings with exact 
solutions of the form X n = s(p m 6), with s(x) being a meromorphic peri- 
odic function. Furthermore, if we try to construct ergodic dynamical systems 
with exact solutions without using addition theorems of elliptic functions, we 
must avoid some lack of analyticity on a map / due to the Weierstrass the- 
orem on functions possessing an algebraic addition theorem as follows. The 
Weierstrass theorem says that any meromorphic function s(u) possessing an 
algebraic addition theorem is either an elliptic function or is of the form R(u) 
or R(e Xu ), where R is a rational function [|l^] . Thus, all of meromorphic and 
periodic functions s(u) possessing an algebraic addition theorem belong to 
the class of elliptic functions. This is the reason why all of our ergodic ra- 
tional transformations constructed here does not have a density function in 
a form p(x) = , 1 = with m > 5 whose integral defines the 

^b +b 1 x+b 2 x 2 + -+b m x m 

inverse function of a hyperelliptic function with genus g(> 2), but has a den- 
sity function in a form p(x) = —, = , whose integral defines 

y b +bix+b2X 2 +b 3 x 3 +b4X i 

the inverse function of an elliptic function with genus g(= 1). 

Of course, not all of ergodic transformations on I are in this category as 
follows: 

Example 1: Gauss Map (1845) 

The Gauss map G : x — > 1/x — [1/x] which maps I — {0} onto I is known 
to be ergodic with respect to the absolutely continuous invariant measure 
li(dx) = t — — — . However, the transformation itself is not a rational function. 

^ v ' lH2(l+a:) ' 

Example 2: Bool Transformation (1857) 

The transformation of Bool B : x — > x — - which maps R — {0} onto R 
is known to be ergodic | 13[] . The transformation preserves infinite measures. 
However, the exact solutions X n do not seem to exist. 
Example 3: 

The two onto mappings S2 '■ x — > 2x mod 1 and S3 : x — >• 3x mod 1 
are mutually permutable in the sense that S3 o S 2 (£) = S2 o Ss(x). They are 
ergodic with respect to the Lebesgue measure on i". However, the transfor- 
mations are not rational. 
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Example 4: 

The two transformations F 2 : x — > \{x — 1/x) and F 3 : x — > 3^231 °f ^ ne rea l 
line i? are mutually permutable in the sense that F3 o F 2 (x) = F2 o ^(x). 
They are ergodic with respect to a unique invariant measure fx(dx) = n ^ x 2^ 
(Cauchy distribution) which is absolutely continuous with respect to Lebesgue 
measure [14[] . These properties comes from the fact that the transformations 
are equivalent to the multiplication formula of cot(x). Thus, these mappings 
have exact solutions X n = cot(2 n #) and X n = cot(3 ra #). Hence, maps of ex- 
ample 4, although they are defined on the infinite support, can be considered 
to belong to the present category of exactly solvable chaos. 



5 Classification problem of exactly solvable 
chaos 

It is known that any elliptic function s(x) can be expressed in terms of Weier- 
strassian elliptic functions p(x) and p'(x) with the same periods, the expres- 
sion being rational in p{x) and linear in g/(a;)|J]. The invariants ^2(^15^2) 
and g$(uJ\, uj 2 ) of Weierstrassian elliptic functions p(x) are not changed under 
a transformation 

u)[ = cnoi + bu 2 , oj' 2 = OjJ\ + du 2 , (52) 

where a,b,c and d are integers satisfying ad — be — 1. An elliptic modular 
function J of p(x) is given by 

jw = A (53) 

where r = Since the elliptic modular function J(r) in Eq.(|53D is invariant 
under a certain condition; namely, 

c^rfr ad -be = I, (54) 

a + or 

where a, 6, c and ci are integers, J(r) of an elliptic function s(x) can serve as 
a characteristic of exactly solvable chaos. For an example, let us consider the 
classification problem of the two-parameter family of the generalized Ulam- 
von Neumann maps X n+1 = /^(X^). In this case, it is easy to check that 

ju m) - m-m)(i-i) + (i- m y? 

J{l,m)- 27(1 _ /)2(1 _ m) 2 (/ _ m) 2 • ^ 

Thus, the equality J(l, m) = J{1\ m') gives us the solutions of the classifica- 
tion problem as follows: 

A = ^A = l-^A = r ^ 7 ,A = i,A = x ^ T ,A = ^, (56) 
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where A = ^— ^ and A' = \ — ^y. This means that when one of the conditions in 
Eq.(|5(]) is satisfied, dynamical systems X n+ \ = //^(X n ) and Y n+ \ = fi? m i{Yn) 

(2) (2) 

have an algebraic relation such that the relation H[fiJ l (X),fi,^ n ,(X)\ = 
holds, where if is a polynomial function in two variables. In this way, we can 
solve the classification problem of exactly solvable chaos X n+1 = f(X n ). 



6 Summary and Discussions 

Rational ergodic maps on the unit interval are systematically constructed 
from addition theorems of elliptic functions in a unified manner. According 
to the classical Weierstrass's theorem about meromorphic functions possess- 
ing algebraic addition theorems and Ritt's theorem about permutable ratio- 
nal functions, such rational mappings due to addition theorems of elliptic 
functions are essentially all rational ergodic transformations that have ex- 
act solutions as well as exact density functions. Therefore, in discrete-time 
systems, solvable models showing chaotic behavior do not belong to an excep- 
tional class but they are ubiquitous as the presence of multiplication theorems 
of elliptic functions. This implies the gap between the notion of solvability 
and integrability in discrete-time dynamical systems as follows: The notion 
of solvability in discrete-time systems can be compatible with chaotic behav- 
ior which are commonly believed to indicate non-integrability of the systems, 
whereas the notion of the solvability is essentially same as that of the integra- 
bility in continuous-time systems (ODE). Thus, further studies of asking what 
is an essential characteristic distinguishing chaotic discrete-time systems from 
regular discrete-time systems are clearly needed for the settlement of the issue. 
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